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ABSTRACT 


The  model  of  a  Surface  Effect  Ship  was  refined  to  include 
simplified  propulsion  dynamics,  negative  drag  characteris¬ 
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sw 


ss 


sp 


Average  sidewall  wetted  area, 
starboard  side 

Average  sidewall  wetted  area, 
port  side 

Added  mass  coefficient  in  roll 
force  equation 

Added  mass  coefficient  in  pitch 
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Average  wetted  sidewall  area 
of  the  bow 
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rad/sec 
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engine 
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engine 
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I.  INTRODUCTION 


Rising  fuel  costs,  critical  manpower  shortages  and  an 

increasing  need  for  rapid  deployment  of  forces  in  defense 

of  United  States'  interests  throughout  the  world  mandate 

the  development  of  fuel  efficient,  high  speed,  minimum  man- 

bb 

ned  ships  for  the  Navy.  The  Captured  Air  Bundle  (CAB)  Sur¬ 
face  Effect  Ship  (SES)  holds  great  promise  in  providing  an 
answer  to  these  problems. 

A  ship  with  minimum  manning  traveling  at  the  high  speeds 
envisioned  for  the  SES  would  require  highly  skilled  operator 
personnel.  It  would  be  necessary  to  provide  these  personnel 
with  real  time  operational  training  similar  to  that  which 
pilots  receive  in  aircraft  simulators.  This  would  test  the 
man  in  responding  to  failure  situations  that  may  be  encoun¬ 
tered  during  actual  operation  and  train  personnel  in  control 
techniques  that  are  unique  to  SES  type  vehicles.  A  Real  Time 
Simulator  could  also  be  used  as  a  means  to  test  hardware  used 
in  SES  craft  and  provide  a  design  tool  for  the  development 
of  future  modifications. 

A  simplified,  non-linear  five  degree  of  freedom  Real  Time 
Simulator  (RTS5D)  for  the  3K-SES  has  been  developed  by  T.  S. 
Nelson  [Ref.  1]  such  that  a  continuously  observable  real  time 
solution  is  interfaced  with  man-generated  control.  Certain 
"man-in-the-loop"  experiments  have  been  conducted  using  the 
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RTS5D  that  have  demonstrated  limitations  in  that  model.  It 
was  necessary  to:  1)  refine  ship  dynamics;  2)  introduce  sim¬ 
plified  propulsion  dynamics;  3)  constrain  thrust  angle  in 
order  to  remain  within  the  bounds  of  programmed  dynamic  res¬ 
ponse  (avoid  broaching) ;  4)  refine  the  software  to  provide 
faster  computer  iteration  time;  and  5)  change  the  hardware 
to  provide  the  operator  with  speed  control  as  well  as  turning 
control . 

Accordingly,  changes  have  been  made  in  the  RTS5D  model 
to  more  accurately  reflect  the  nonlinear  forward  drag 
characteristics  of  the  SES.  A  speed  controller  has  been 
designed  to  allow  operation  of  the  3K-SES  within  a  speed 
range  of  40-60  knots  and  limits  have  been  placed  on  the 
thruster  deflection  angle  in  order  to  prevent  broaching,  a 
condition  which  is  not  allowed  under  the  existing  dynamic 
model.  Software  refinements  have  been  incorporated  to  im¬ 
prove  computer  iteration  time. 


II.  EQUATIONS  OF  MOTION 


A.  COORDINATE  SYSTEMS  AND  ASSUMPTIONS 

The  RTS5D  model  of  the  3K-SES  used  simplified  equations  of 
motion.  For  completeness  of  comparison  the  development  of 
these  equations  is  repeated  here.  The  rational  for  requiring 
simplification  was  to  start  with  differential  equations  of 
motion  that  could  be  utilized  effectively  in  a  real  time  simu¬ 
lation.  To  satisfy  this  requirement  a  model  was  selected  that 
would  use  point  source  of  drag  forces  consisting  of  only  a 
few  lumped-coefficient  terms.  The  equations  used  are  based 
upon  those  developed  in  Ref.  2  which  follow  this  concept. 
Additional  assumptions  used  in  the  simplified  5  DOF  equations 
were  as  follows: 

1)  All  accelerations  are  measured  at  the  center  of 
gravity. 

2)  Vehicle  is  "free- to-heave" . 

3)  All  cross  coupled  moments  of  inertia  are  zero. 

4)  All  cross  products  of  angular  velocities  in  force 
equations  are  zero. 

5)  All  roll  and  pitch  angles  used  in  force  equations 
are  subject  to  small  angle  approximations. 

6)  Calm  water  conditions. 

7)  Mass  and  mass  distribution  of  vehicle  are  constant. 

8)  Effect  of  changes  of  aerodynamic  force  are  neglected. 
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Utilizing  the  above  assumptions,  the  simplified  S  DOF 
equations  of  motion  are  developed  using  the  coordinate  sys¬ 
tems  shown  in  Fig.  1  and  Fig.  2,  the  force  vector  diagrams 
shown  in  Fig.  3,  Fig.  4A,  Fig.  4B,  and  the  following 
equations  from  Newton’s  laws  of  motion: 


SURGE 

m(u  - 

VT)  -  X 

SWAY 

m(v  - 

ur)  -  Y 

YAW 

1  z  y 

*  N 

PITCH 

I  y  q 

3  M 

ROLL 

rx  p 

*  K 

where : 

m 

mass  of 

the  rigid  ship 

'  slugs 

u  * 

velocity 

in  x-direction  (SURGE) 

'  ft/sec 

v  =  velocity  in  y-direction  (SWAY)  ’  ft/sec 

r  a  angular  velocity  about  the  z-axis  '  deg/sec 

p  =  angular  velocity  about  the  x-axis  '  deg/sec 

q  =  angular  velocity  abouc  the  y-axis  '  deg/sec 

? 

I-  =  moment  of  inertia  about  the  z-axis  '  lb  -ft 

l  m 

7 

Iv  =  moment  of  inertia  about  the  y-axis  *  lb  -ft 

x  m 

2 

Iv  =  moment  of  inertia  about  the  x-axis  ’  lb  -ft 

a  m 

X  =  summation  of  forces  in  x-direction  '  lb^ 

Y  *  summation  of  forces  in  y-direction  '  lb^ 

N  =  summation  of  moments  about  the  ’  ft-lb- 

z-axis  * 

M  »  summation  of  moments  about  the  y-axis  ’  ft-lb^ 

K  *  summation  of  moments  about  the  x-axis  '  ft-lb^ 


Additionally  the  following  navigation  relationships  are 


utilized  (see  Fig.  1). 


X 

0 

■  U  CO  s  ^ 

-  v  sinn» 

*0 

■  u  sin^ 

♦  V  COSij; 

V  2 

s 

-  u2  ♦  V2 

8 

■  tan  * 

(-v/uj 

where : 

V  *  heading  angle 
8  -  drift  angle 
5  -  thrust  vector  angle 

Also , 


X 

0 

* 

n  dt 

0 

u 

n 

/  U 

dt 

Y 

0 

* 

/Y  dt 

0 

V 

- 

/V 

dt 

r 

= 

/r  dt 

♦ 

- 

;r 

dt 

P 

- 

/p  dt 

* 

X 

fP 

dt 

q 

m 

/q  dt 

0 

X 

f  q 

dt 
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N'ote  : 


A  :;a v 


ynav 


1)  ^=0  when  vehicle  heads  paralled  to  ♦  XV4V 

axis.  NAV 

2)  Vehicle  in  figure  above  shown  with  negative 
thrust  vector  angle  6,  positive  ,  negative 
6,  negative  sway  velocity  v,  i.e.,  in  a  right 
turn 


Figure  1 

Definition  Of  Coordinate  System  (Part  I) 
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Note:  Direction  for  positive  velocity  and 
angular  rates  are  shown 

Figure  2 

Definition  Of  Coordinate  System  (Part  II) 
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1)  p  and  p'  are  equivalent  point  force  centroids 

2)  <5^ ,  5 g ,  are  thrust  vector  angles  and 

are  not  required  to  be  equal. 

3)  T^,  Tg,  Tg ,  Tjq  are  thrust  magnitudes  and  are 

not  required  to  be  equal. 

4)  Force  directions  shown  are  for  a  right  turn 


Figure  3 

Surface  Effect  Ship  (Top  View) 


B.  FORCES  AND  MOMENTS 

1.  Surge  Forces  (see  Fig-  4B) 

The  forward  acceleration  was  determined  to  be  genera¬ 
ted  by  the  summation  of  thrust  vectored  in  the  forward  direc¬ 
tion  T£o  .  The  resultant  acceleration  of  the  SES  was 
assumed  to  be  opposed  by  a  retarding  force  exhibiting  a 
"velocity- squared"  characteristic.  Additionally,  Newton's 
laws  of  motion  specified  a  "v.r"  product  that  contributed 
to  the  drag  component  of  the  equation  when  the  vehicle  was  in 
a  turn  (contrifugal  force  reaction  term).  The  simplified 
equation  of  motion  for  the  surge  acceleration  follows  from 
a  summation  of  forces  in  the  u  direction  as  depicted  in 
Fig.  4B. 

u  =  (Tforw/m)  ‘  (CDXu2/m)  +  vr 

2.  Sway  Forces  (see  Fig.  4A) 

The  side  acceleration  of  the  vehicle  was  analyzed 
to  be  the  result  of  the  summation  of  the  thrust  vectored 
parallel  to  the  vehicle's  y-axis  and  a  retardation  "velocity 
squared"  phenomena  such  as  that  produced  by  the  cross  flow 
drag  term  of  the  sidewall  as  described  in  Ref.  3.  This 
retardation  force  is  augmented  by  a  "u-r"  product  as  speci¬ 
fied  by  Newton's  Law. 

v  =  Tside/m  •  CDY  v|v|/m  '  ur 
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Note  the  requirement  to  force  the  C^yV^/m  term  to 
maintain  the  sign  of  the  sway  velocity  such  that  both  left 
and  right  turns  may  be  computed. 

3.  Yaw  Moments  (see  Fig.  3) 

The  yaw  acceleration  was  a  turning  moment  summation 
which  was  defined  to  be  influenced  by  three  primary  forces 
operating  over  moment  arms.  It  is  of  interest  to  note  that 
the  third  term  is  an  "added  mass"  term  which  is  considered 
to  operate  on  the  same  pressure  point  as  the  cDyvlvl  term. 
The  term  T^aw  is  a  summation  of  moments  generated  by  the 
four  thrusters  as  shown  in  Fig.  3. 


Tyaw  ~  S4T7  cos<S7  +  s3T8  cos  58  "  S1T9  C0S  59 


-  s~Tir.  cos  6 .  n  -  l  T_  sin  5_  -  l  T0  sin  S0 
Z  1 0  1U  O  7  7  Oo  o 


*oT9  Sin  S9  "  2'oT10  sin  510 


r  =  CDY  4w  VIVI/IZ  +  A22UVVX2  +  TYAW/XZ 

4.  Pitch  Moments  (see  Fig.  4B) 

The  pitch  acceleration  was  assumed  to  be  the  summa¬ 
tion  of  moments  generated  by  forward  thrust,  f£orw,  the 
buoyancy  of  the  sidewalls  of  the  SES,  F^  and  F^  and  a 
vertical  force  generated  at  the  bow  of  the  vehicle,  F^. 

This  vertical  bow  force  was  defined  to  be  a  lumped  parameter 
term  which  modeled  the  reaction  force  due  to  plenum  pressure 
acting  against  the  bow  seal. 
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=  average  draft  of  bow  seal 


A 


w. 


Vx  +  ((ix  tan9^l/2)1x 


A 


w- 


A 


w, 


X2 


'  Vxj  ■  (fix  tane)/nx 


=  A  5,  p  g 
W1  X2 


'xl 


Aw 9  !'x, 

CDpPbWel'td  *  *31  tan  9") 

=  average  draft  of  sidewall 
=  plenum  pressure 

=  average  wetted  sidewall  area  of  the  stern 
=  width  of  one  sidewall 


w 


w- 


A33 


p 

g 


'31 

'3 


width  of  bow  seal 

average  wetted  sidewall  area  of  the  bow 

added  mass  coefficient 
water  density 
32.3  ft/s2 

lever  arm  of  bow  seal 
lever  arm  of  buoyancy  force 


q  =  * (Tr  l  +  +  F-l-  -  C_  u“  i 

n  v  forw  p  p  31  2  p  DX  ; 


Fl*3  *  A33Uq')/IY 
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Note  the  added  mass  term  A^uq  which  was  required 
to  provide  damping  to  the  pitch  moment  and  is  specified 
in  Ref,  2. 


It  was  found  that  the  response  of  the  SES  to  a 

pitch  perturbation  without  the  added  mass  component  was 

undamped  and  approximately  sinusoidal.  Unlike  the  yaw 

equation  where  the  added  mass  term  was  a  small  contributor 

to  damping,  the  pitch  added  mass  term  was  found  to  be  of 

significant  importance  in  modeling  the  known  pitch  motion. 

S.  Roll  Moments  (see  Fig.  4A) 

The  roll  acceleration  equation  was  analyzed  to  be  a 

summation  of  moments  generated  by  buoyancy  forces  of  the 

port  and  starboard  sidewalls,  F  and  F  ,  the  thrust 

sp  ss 

vectored  parallel  to  the  y-axis  of  the  vehicle,  T  a  side 

force  v|vj,  and  a  lumped  coefficient  vertical  force, 

F  .  The  vertical  force  F  was  defined  as  the  force 
sw  sw 

generated  in  a  turn  due  to  the  sidewall  curvature  (dead 
rise  angle)  acting  against  the  cross  flow  of  water. 
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A^p  =  average  wetted  area  port 
A^s  =  average  wetted  area  starboard 


F  )l 
ss  sw 


T  i  + 
side  p 


Z 

z 


-  F  l  -  A,,  up)/I_, 
sw  sw  31  VJ  X 

Again  note  the  added  mass  term  A^up/I^  which  was 
found  to  be  essential  in  modeling  the  damping  phenomena  of 
the  roll  motion. 

In  summary,  the  simplified  5  DOF  equations  of  motion 


used  in 

RTS5D 

for  the  3K  TON  SES 

are : 

SURGE 

u  = 

Tforw/ra  -  CDX  u2/m 

+  vr 

(Fig. 

4B) 

SWAY 

v  = 

Tside/ra  *  CDY  V'V| 

/m  -  ur 

(Fig. 

4A) 

YAW 

r  - 

Tyav/XZ  +  CDY 

V1  z 

*  A22uvtw/Iz 

(Fig. 

3) 

PITCH 

q  = 

(Tforw*p  +  F3*31  + 

F2*3 

7 

Cpy  u“lz 

-  F1a3  -  A34uq)/Iy  (Fig.  4B) 

P  3  CCFsp  -  Fss).sw  -  Tside*p  ♦  CDy  v|vU2 

'  Fsw*sw  -  A33  up)/IX  (Fig-  4A) 
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C.  PARAMETER  IDENTIFICATION 

These  equations  are  an  extension  of  the  3  DOF  flat  turn 
SE5  model  developed  by  Gerba  and  Thaler  in  Ref.  4.  The 
identification  of  craft  parameters  C^,  C^y,  and  i  is 
described  in  Ref.  4  and  repeated  here  for  completeness. 

The  surge  drag  coefficient  is  determined  by  selecting 
a  steady  state  turn  condition,  where  u,  v,  r,  and  m 

are  known. 

i  -  0  -  -  Ss“!  -  vr 


from  which 
CDX 


m 


Tforw  *  ”vr 


m 


The  sway  drag  coefficient  C^y  is  determined  by  the  same 
steady  state  turn  condition  which  requires 
T  C 

n  side  DY  ,  , 
v  *  0  =  -  -  — —  v  v  -  ur 


m 


m 


where 


'DY 


,  Tside  ;  mur 
vlv! 


The  sway  drag  moment  arm  follows  utilizing  the  yaw  accel¬ 
eration  equation 

•  „  -  Tyaw  „  CDYv|v|jw  „  A22  uv2w 

r  =  0  =  -j—  ♦  - T -  +  - 1 - 

L1  Ll  LZ 
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%  — 


i 


which  yields 


w 


-T 

yaw _ 

v|v|  +  A^uv 


The  addition  of  the  pitch  and  roll  equations  introduces 

additional  lumped  coefficients  C^p  and  C^p.  These  are 

easily  solved  using  known  constants  l  ,  1  ,  Cnv,  z  ,  z  , 

1  &  w*  p*  DY*  z  sw* 

w  ,  Pw>  1 a  steady  state  values  of  F  ,  F  ,  F.,  F,, 

e’rbd  /  spss’1’2’ 

Tf ,  T  .  .  ,  u,  v,  and  r.  It  is  significant  to  note  that 

in  a  steady  state  condition  the  angular  roll  rate,  p,  and 
angular  pitch  rate,  q,  are  both  required  to  equal  zero. 

This  is  in  contrast  to  the  yaw  rate,  r,  where  a  finite 
steady  state  value  is  desired  in  a  turn.  Thus  the  added 
mass  terms  of  A^-uq  and  up  in  the  pitch  and  roll  equa¬ 
tions  are  not  utilized  in  the  determination  of  craft  para¬ 
meters;  their  function  is  strictly  confined  to  damping  of 
their  respective  accelerations.  Therefore,  utilizing  the 
pitch  equation  with  known  steady  state  turn  values. 


q  =  0  tTforw*p  +  CDPpbWe(-Jld  '  *31  tan  0^*31 
+  F2*3  *  CDXU  *2  ’  F1*3^IY 


yields 

C 


DP 


CDXU  * z  *  Fl*3  ~  F2*3  "  Tforw  *p 
pbV*d  *  *31  tan  0)*31 
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III.  RTS5D  VALIDATION 


A.  DBSIM5D  BENCHMARK 

Validation  of  the  RTS5D  was  attempted  using  the  data 
base  program  (DBSIM5D)  [Ref.  5]  as  a  benchmark.  However, 
new  guidance  in  the  proper  use  of  the  data  base  program  con¬ 
tained  in  Ref.  6  revealed  the  following  errors  in  the  use  of 
the  DBS IMS D.  1)  Initial  conditions  had  been  improperly  set 
and  2)  Data  used  for  validation  for  the  56  knot  run  at 
effector  deflection  angles  greater  than  10°  was  not  correct 
because  of  the  thrust  inlet  broach  conditions  which  existed 
using  the  DBSIM5D.  These  conditions  were  not  accounted  for 
in  the  RTS5D  model.  The  noted  discrepancies  in  the  opera¬ 
tion  of  the  DBSIM5D  necessitated  a  revalidation  of  the 
RTSSD. 

B.  VALIDATION  RESULTS 

,  Using  the  information  contained  in  Ref.  6,  the  DBSIM5D 
program  was  again  used  as  a  benchmark  to  validate  the  RTSSD 
simulation.  A  DBSIM5D  sequence  of  runs  with  initial  forward 
velocities  equal  to  40  knots,  50  knots  and  60  knots  with 
step  effector  angles  of  5°  and  10°  (at  40  knots),  5°,  10° 
and  15°  (at  50  knots) ,  and  5°  (at  60  knots)  were  compared 
to  identical  runs  of  the  RTSSD.  Thrust  effector  angles 
were  kept  below  a  maximum  that  would  have  caused  a  broach 
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condition.  Broach  conditions  are  explained  in  Chapter  IV. 

The  first  peak  overshoot  value  and  "Quasi  Steady  State" 
values  for  the  variables  u  (surge),  v  (sway),  r  (yaw  rate), 

<j>  (roll),  and  0  (pitch)  were  used  as  a  basis  for  comparison. 
The  results  are  shown  in  Tables  I  - 1 1 1 .  Significant  error 
in  forward  velocity  occurred  in  the  50  knot  test.  Also,  it 
was  noted  that  in  40  knot  test  with  5°  thruster  deflection 
angle  that  the  DBSIM5D  actually  reached  a  "Quasi  Steady 
State"  forward  velocity  greater  than  the  initial  value  upon 
entering  the  turn.  These  conditions  led  to  the  development 
of  the  model  of  the  RTSSD  with  negative  drag  characteristics 
in  sea  state.  Differences  in  the  other  measured  parameters 
were  addressed  by  parameter  adjustment  on  the  new  model  in 
ordeT  to  more  closely  match  values  obtained  with  the  DBSIM5D. 


TABLE  I 


RTS5D  and  DBSIM5D  Performance  Test  at  40  Knots 

effector 

angle  DBSIM5D  RTS5D 

<j>  in  degrees  <j> 


1st  pk 

Qss 

1st  pk 

Qss 

5° 

1.01 

.  742 

.17 

.12 

10° 

1.46 

1.39 

.31 

.24 

1st  pk 

* 

in 

Qss 

degrees 

1st  pk 

* 

Qss 

5° 

1.44 

1.40 

1.27 

1.21 

10° 

1.48 

1.48 

1.28 

1.20 

1st  pk 

u 

in 

Qss 

ft/sec 

1st  pk 

u 

Qss 

5° 

67.4 

70.9 

n/a 

67.03 

10° 

n/a 

66.3 

n/a 

65.98 

1st  pk 

V 

in 

Qss 

ft/sec 

1st  pk 

V 

Qss 

5° 

1.42 

1.37 

2.84 

2.42 

10° 

3.20 

3.26 

3.83 

3.41 

1st  pk 

r 

in 

Qss 

degrees/sec 
1st  pk 

r 

Qss 

5° 

.396 

.403 

.  58 

.33 

10° 

.953 

.959 

.98 

.  66 
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TABLE  II 


RTS5D  and  DBSIM5D 

Performance 

Test  at  50 

Knots 

effector 

angle 

1st  pk 

DBSIM5D 

* 

in  degrees 
Qss 

RTS5D 

* 

1st  pk 

Qss 

5° 

.424 

.446 

.  27 

.19 

10° 

.757 

.  886 

.49 

.38 

15° 

1.02 

1.4 

.72 

.  57 

1st  pk 

e 

in  degrees 
Qss 

9 

1st  pk 

Qss 

5° 

.  862 

.892 

1.16 

1.15 

O 

o 

r—i 

.943 

1.13 

1.15 

1.16 

15° 

1.12 

1.35 

1.22 

1.17 

1st  pk 

u 

in  ft/sec 
Qss 

u 

1st  pk 

Qss 

5° 

n/a 

83.6 

n/ a 

83.8 

10° 

n/a 

78.7 

n/a 

32.51 

15° 

n/a 

69.2 

n/a 

80.58 

1st  pk 

V 

in  ft/sec 
Qss 

V 

1st  pk 

Qss 

5° 

1.31 

1.30 

3.55 

3.02 

10° 

2.72 

2.70 

4.79 

4.26 

15° 

4.05 

4.33 

5.69 

5.207 

1st  pk 

r 

in  degrees/sec  r 

Qss  1st  pk 

Qss 

5° 

.368 

.365 

.  73 

.41 

10° 

.  864 

.870 

1.22 

.  82 

15° 

1.35 

1.44 

1.66 

1.26 
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TABLE  III 


RTS5D  and  SBSIM5D  Performance  Test  at  50  Knots 

effector 


angle 

DBSIM5D 

RTS5D 

in  degrees 

r 

1st  pk 

Qss 

1st  pk 

Qss 

.141 

.200 

.39 

.274 

0 

in  degrees 

0 

1st  pk 

Qss 

1st  pk 

Qss 

5° 

.  321 

.  583 

1.14 

1.07 

u 

in  ft/sec 

u 

1st  pk 

Qss 

1st  pk 

Qss 

5° 

n/a 

97.8 

n/a 

100.29 

V 

in  ft/sec 

V 

1st  pk 

Qss 

1st  pk 

Qss 

5° 

2.03 

1.84 

4.26 

3.62 

r 

in  degrees/sec 

r 

1st  pk 

Qss 

1st  pk 

Qss 

5° 

.532 

.  494 

.  87 

.487 
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IV.  MODEL  DEVELOPMENT 


A.  DRAG  CHARACTERISTICS 

1 .  Surge  Equation  Model 

The  RTS5D  model  assumed  a  velocity  squared  drag  re¬ 
tarding  characteristic  for  the  3K-SES  which  led  to  the  equa¬ 
tion  for  surge  acceleration  that  was  presented  in  Chapter  II. 
The  drag  characteristics  from  Ref.  7  for  the  3K-SES  that  is 
simulated  by  the  DBSIM5D  are  shown  in  Fig.  5.  The  data  base 
program  was  run  using  parameters  that  caused  it  to  exhibit 
a  retarding  force  corresponding  to  the  curve  at  sea  state 
four,  therefore  this  curve  was  used  to  remodel  the  RTS5D. 

The  new  simplified  equation  of  motion  for  surge  acceleration 
became 


(Tforw/nn  '  +  v  •  r 


where  FSD  is  full  scale  drag  which  is  defined  by  the  curve 
in  Fig.  5.  Fpqrw  *s  the  sumation  oF  the  four  thrusts  assumed 
in  this  model  to  be  applied  at  the  same  point  along  the 
center  line  of  the  craft. 

2 .  Implementation 

FSD  was  implemented  into  the  RTS5D  model  by  dividing 


the  curve  into  four  separate  regions  and  defining  FSD  within 
each  region.  In  region  I  the  retarding  force  was  approxi¬ 
mated  by  a  normalized  velocity  cubed  term.  In  regions  II 


VT\ 

I 

o 
** 
ctf  X 
U  W 
Q  & 


r> 


CVi 
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V 


and  III  straight  line  approximations  were  used.  In  region 
IV  a  polynomial  curve  fit  was  used.  The  equations  defined  for 
each  region  are  listed  below. 


Region  I  -  (  0-18  knots) 

PSD  =  CDX1  (  u/  umax  )3 
Region  II  -  (18-22  knots) 

FSD  =  CONSTANT^  -  CDX5  u 
Region  III  -  (22-36  knots) 

FSD  =  CDX2  u  -  CONSTANT^ 
Region  IV  -  (above  36  knots) 

FSD  =  CDX2  u2  +  CDXi+/  u1,3 


3 .  Identification  of  Surge  Drag  Coefficients 

In  region  I,  is  equal  to  the  first  peak  value 

of  drag.  In  regions  II  and  III,  and  Cj^  are  equal  to 

the  slopes  of  the  drag  curve  and  the  constants  were  solved 
for  by  applying  the  boundary  conditions  at  the  intersections 
of  the  adjacent  regions.  The  values  of  CpX,  and  were 

computed  by  taking  FSD  at  two  different  known  values  of  u 
and  solving  for  and  simultaneously. 


B.  LINEARIZED  MODEL  ANALYSIS 


Since  the  surge  model  of  the  RTS5D  is  extremely  non¬ 
linear  over  the  complete  speed  range,  analysis  is  very  com¬ 
plex.  Therefore  the  model  was  linearized  and  it’s  performance 


analyzed  at  various  points  on  the  FSD  curve.  Linearization 
of  the  model  was  accomplished  by  computing  the  linear 
coefficients  through  partial  differention  with  respect  to 
surge  in  each  of  the  four  region  equations  and  evaluating 
these  terms  at  various  operating  point  surge  velocities  (uq) . 
The  linearized  model  is  shown  in  Fig.  6  and  the  linear  coef¬ 
ficients  for  the  surge  drag  force  FSD  are  shown  below. 

3  CDX1  <  V  W2 
"  CDX5 
°DX2 

2  CDX3  Uo  "  1,5  CDX^/  u02’5 

A  map  of  the  linearized  surge  directional  roots  are 
shown  in  Fig.  7.  The  characteristic  roots  for  each  operating 
region  are  shown  coded  (see  legend) .  Note  that  only  Region  I 
and  Region  IV  roots  are  uq  dependent  and  therefore  move  as  uq 
changes.  The  direction  of  change  along  the  real  axis  is 
noted  by  the  arrows  shown.  Region  II  and  III  roots  are  fixed 
in  locations  and  independent  of  uq . 


Region 


/ 


I 

FSD(u0) 

II 

FSD(uo) 

III 

FsD(uo) 

IV 

FSD(uo) 
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Region  IV  of  Fig.  7  is  of  particular  interest  since  the 
ability  to  operate  in  this  speed  range  is  a  primary  consi¬ 
deration  in  the  development  of  the  SES.  It  can  be  seen 
from  the  map  of  the  roots  that  within  this  region,  at  speeds 
from  36  to  almost  60  knots,  poles  exist  in  the  RHP.  This 
is  a  result  of  the  negative  drag  slope  where  the  retarding 
force  actually  decreases  with  an  increase  in  speed,  the 
magnitude  of  variation  being  dependent  on  sea  state.  Region 
II  also  exhibits  this  negative  drag  characteristic  but  since 
it  is  of  such  a  narrow  speed  range,  operation  in  this  re¬ 
gion  was  not  considered  in  this  design.  Further  modeling 
changes  for  the  other  degrees  of  freedom  were  based  on 
Region  IV  response.  In  Regions  I,  III  and  in  Region  IV  above 
60  knots,  the  system  will  be  stable  at  the  speed  uq  where  FSD 
is  equal  to  commanded  thrust. 
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Figure  7 

Linearized  Surge  Root  Map 


1 


If  it  were  necessary  to  operate  at  a  speed  within  the 

negative  drag  slope  characteristic,  it  would  be  necessary  to 

actively  modulate  the  thrust  command  in  order  to  maintain 

* 

a  nearly  constant  speed,  but  this  would  require  a  great  deal 
of  operator  attention  since  the  negative  drag  slope  in  Region 
IV  is  one  of  slowly  divergerrt  instability.  For  this  reason 
a  closed  loop  speed  controller  has  been  included  in  the  data 
base  model  and  therefore  was  also  added  to  the  RTS5D  model. 

C.  SIMPLIFIED  PROPULSION  DYNAMICS 

Before  adding  a  speed  autopilot  to  the  model  one  addi¬ 
tional  refinement  was  included  in  the  design.  Thrust  command 
changes  in  the  RTS5D  as  reported  in  Ref.  I  were  instantan¬ 
eously  entered  into  che  equations  of  motion  as  thrust 
changes.  This  was  not  the  case  for  the  data  base  simulation 
program.  Specific  thrust  levels  are  generated  by  specific 
water  flow  rates  through  water  jet  nozzles.  These  flow 
rates  are  nearly  proportional  t-o  the  pump  speed  that  pro¬ 
duces  them.  Since  the  pumps  are  driven  directly  by  gas 
turbine  engines,  a  change  in  thrust  occurs  with  a  change  in 
turbine  speed.  A  thrust  command  translates  to  a  power  tur¬ 
bine  speed  command  which  as  an  output  has  power  turbine 
speed  achieved.  This  output  is  translated  into  thrust 
achieved.  It  was  specified  in  Ref.  7  that  one  of  the  goals 
in  the  base  line  design  was  to  have  power  turbine  speed 
control  response  c harast erist ics  which  are  relatively  inde¬ 
pendent  of  the  power  setting.  Accordingly,  in  the  design 
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of  the  auto  pilot  of  Ref.  7,  the  simplified  propuLsion 
block  shown  in  Fig.  8  was  used.  The  magnitude  of  the  pro 
pulsion  system  time  constant,  t  was  not  specified.  For 
the  purpose  of  this  analysis  r  was  chosen  to  be  1  second. 


Command 

1 

IT 

l  +Ts 

Figure  8 

Simplified  Propulsion  System 


D.  ANALYSIS  OF  LINEARIZED  MODEL  WITH  SPEED  CONTROLLER 

The  complete  linearized  model  with  speed  controller  that 
is  included  in  the  model  for  the  purposes  of  this  study  is 
shown  in  Fig.  9.  The  saturation,  effects  of  fuel  control  are 
ignored  and  the  linear  representation  is  shown  as  which 
is  assumed  to  be  unity.  The  controller  includes  propor¬ 
tional  plus  integral  control.  The  proportional  control  was 
included  to  keep  the  ship  speed  error  low  while  t.he  integral 
control  was  included  in  order  to  tero  out  long  term  control 
error . 
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Automatic  Thrust 


Propulsion  System  and  Speed  Controller 


The  characteristic  equation  for  the  system  in  Fig.  9  is 


S3  +  (1  +  FSD/  J  m  )  S2  +  l/m  (PSD/  ,  +  K)  S  +  K/  mf  s° 
vuo'  '  o' 

Using  the  Routh  Hurwitz  Criterion  for  stability 
s3  1  1/m  <FSD#(  }  +  K) 

S2  (  1  +  FSD1/  ./  m  )  K/  m  T 

v  o; 

S1  1/m  (1  +  FSD'(u  j/m)  (FSd'(u  j  +  K)  -  K /bT  0 

s0  K/  m  r  0 

For  Stability  K  must  be  greater  than  0  and  greater  than 

PSD'(u0)  T  +  ”rFSD'(u0) 

^(uj^  +  —  m 

o 

The  values  of  controller  parameters  K  and  t  can  be  initialized 
by  the  operator  for  either  loosely  controlled  speed  or  tightly 
controlled  speed  for  such  situations  as  station  keeping, 

UNREPS,  or  test  operations.  For  comparison  to  the  data  base, 
these  parameters  were  selected  to  achieve  closest  agreement 
between  the  values  of  surge  achieved  by  the  RTS5D  and  DBSIMSD 
models  after  the  completion  of  a  360°  turn  with  50  knot  initial 
surge  veloctiy  and  a  thrust  effector  angle  of  IS8.  This  test 


condition  for  the  nonlinear  models  was  chosen  for  two  reasons 
First,  given  fuel  conservation  as  a  necessity,  50  knots  is  a 
likely  high  speed  initial  condition  because  it  is  close  to 
minimum  drag.  Second,  this  initial  condition  provided  an 
opportunity  to  test  the  operation  of  the  speed  controller  in 
the  negative  drag  region  of  FSD  since  the  decrease  in  speed 
going  into  the  turn  with  uq=50  knots  caused  increased  drag. 
The  values  of  controller  parameters  that  caused  closest  agree 
ment  between  the  two  nonlinear  models,  when  tested  in  the 
linearized  model  showed  instability  in  the  linear  sense. 

This  is  desirable  in  a  turning  mode  for  the  nonlinear  model 
because  it  allows  speed  to  decrease  and  the  turn  to  be  com¬ 
pleted  more  quickly.  Stiff  speed  control  could  be  achieved 
by  increasing  the  controller  gains,  but  would  normally  be 
used  only  in  straight  ahead  runs. 

Coupling  of  the  equations  of  motion  developed  in  Chapter 
II  necessitated  a  repetition  of  the  parameter  identification 
procedure  resulting  in  new  values  for  C^^ ,  tw,  CDp  and  CD7p. 
These  values  were  adjusted  for  the  best  curve  fit  compared 
to  the  data  base  using  the  linearized  model  for  the  RTS5D 
in  a  360°  turn  using  15°  thrust  effector  angle  with  initial 
surge  velocity  50  knots.  An  identical  series  of  tests  to 
those  in  Chapter  III  were  conducted  comparing  the  new  model 
of  the  RTS5D  to  the  DBSIM5D.  Complete  validation  results 
are  tabulated  in  Chapter  VI. 


E.  EFFECTIVE  THRUST  EFFECTOR  ANGLE 


As  a  result  of  the  negative  drag  characteristics  imposed 
on  the  model,  RTS5D  yaw  rate  (r)  was  high  in  the  low  speed 
test  and  low  in  the  high  speed  test.  This  opposite  error 
at  both  ends  of  the  test  speed  spectrum  was  reduced  by 
introducing  a  correction  factor  to  the  thrust  effector  angle 
proportional  to  yaw  rate 

2  --  =  Z  -  SLIP  •  r 
eff 

and  adjusting  £w  to  increase  r.  This  had  the  desired  effect 
of  reducing  error  at  both  ends  of  the  test  spectrum  since 
the  reduction  of  the  thrust  effector  angle  was  greater  at 
the  high  end. 

F.  BROACH  CONDITION  FLAG 

As  stated  in  Chapter  I,  broach  response  characteristics 
are  not  included  in  this  model,  therefore  a  broach  condition 
warning  flag  was  included  in  the  design  in  order  to  alert 
the  operator  who  could  then  take  corrective  action  to  avoid 
broaching. 

Broaching  is  caused  by  a  water  intake  being  exposed  to 
air  and  is  a  function  of  surge  speed,  roll  angle,  pitcn 
angle,  drift  angle  and  sea  state.  For  the  data  base  3K-SES 
with  broach  boundaries  defined  at  10%  air  by  volume,  the 
broaching  boundaries  in  calm  water  in  terms  of  roll  angle, 
pitch  angle  and  drift  angle  are  shown  in  Fig.  11  for  40  knots 
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and  Fig.  12  for  60  knots.  Reference  6  listed  maximum  thrust 
effector  angles  as  a  function  of  surge  speed  and  air  flow 
in  order  to  avoid  broaching.  This  condition  is  shown  plotted 
in  Fig.  10  for  fixed  plenum  air  flow  rate. 


EFFECTOR 

ANGLE 


Figure  10 

Maximum  Thrust  Effector  Angle  As  A  Function 
Of  Surge  Velocity  In  Order  To  Avoid  Broaching 


As  stated  in  Chapter  III  the  model  developed  in  Ref.  1  did 
not  flag  the  broaching  condition.  To  implement  the  broach 
warning  flag  and  display  the  other  new  parameters  for  the 
modified  RTS5D  the  pilot  graphic  display  described  in  Ref.  1 
was  changed  to  incorporate  the  necessary  data.  Figure  13A 
depicts  typical  information  displayed  to  the  pilot  under 
normal  operating  conditions.  The  new  information  to  the 
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operator  shown  in  the  lower  right  of  Fig.  1 3A  is  sea  state 
(SEAST) ,  computer  assist  (CASST,  on=l  off=0j,  controller 


gains  (K,KK)  and  speed  commanded  (SPCOM) .  If  the  operator 
exceeds  the  thrust  effector  limits  (IMAX)  to  prevent 
broaching  (see  Fig.  10),  a  flashing  warning  is  displayed 
in  place  of  the  bottom  line  of  information  on  the  screen. 
This  is  shown  in  Fig.  13B. 
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Figure  12 

Broaching  Boundaries  at  60  Knots 


VARIABLE  DYNAMICS 

U  (HEAVY) 

UDOT( LIGHT) 


V  (HEAVY) 
VDOT( LIGHT) 


R  (HEAVY) 
RDOT (LIGHT) 


IE 


P  (HEAVY) 
PDOT( LIGHT) 


Q  (HEAVY) 
QDOT (LIGHT) 


yen* 


CONTROL  INPUTS  TIME  NAVIGATIONAL  DATA 

THRUST  RUDDER  DELT  ELAPSED  PITCH  HEADING/RATE  DRIFT  ROLL  SPEED 

300000  -.20  .084  22.2  .52  325/2.1  -.28  .8  50.1 

T7  T8  T9  T10  SEAST  CASST  S  KK  SPCOM  ' 

75000  75000  75000  750000  4  1  2 900  3  50 


Figure  13A 

Pilot  Graphic  Display  (Normal) 


VARIABLE  DYNAMICS 


SEA  TRACK 


U  (HEAVY) 
UDCT( LIGHT) 


V  (HEAVY) 
VDOT( LIGHT) 


R  (HEAVY)  HC 
RDOT( LIGHT) 


P  (HEAVY) 
PDOT( LIGHT) 

Q  (HEAVY) 
QDCT( LIGHT) 


I 


CONTROL  INPUTS  TIME  NAVIGATIONAL  DATA 

THRUST  RUDDER  DELT  ELAPSED  PITCH  HEADING/RATE  DRIFT  ROLL  SPEED 

3000000  -.52  .084  22.2  .52  325/2.1  -.28  .8  50.1 

T7  T8  T9  T10  SEAST  CASST  K  KK  SPCOM 

—  B  R  0  A  C  H  — 

Figure  13B 

Pilot  Graphic  Display  (Broach  Condition  Exists) 


53 


V.  IMPLEMENTATION 


A.  INTRODUCTION 

Given  the  existence  of  a  Real  Time  Simulation  of  a  SES, 
a  decision  was  made  to  refine  that  model  using  more  accurate 
equations  of  motion,  improved  computer  iteration  time  and 
operator  hardware  modifications.  A  complete  description  of 
the  RTS5D  is  contained  in  Ref.  1,  therefore  only  changes 
are  discussed  here,  and  included  are  the  introduction  of 
negative  drag  in  a  sea  state,  propulsion  dynamics,  speed 
control  and  a  broach  condition  warning  in  software.  Hard¬ 
ware  modification  included  changing  the  thruster  control 
box  to  make  it  more  accessable  to  an  operator. 

B.  REQUIREMENTS 

The  following  criteria  were  required  of  the  remodeled 
RTS5D. 

1)  The  refined  simplified  equations  of  motion  would  be 
solved  and  the  results  output  on  a  real  time  basis. 

2)  The  solution  would  be  subject  to  real  time  control 
efforts  generated  by  an  operator  observing  the  output. 

3}  Computer  iteration  time  was  required  to  be  less  t-han 
10  0  ms. 

A  more  complex  discussion  of  condition  3)  is  presented  in 
Appendix  A. 
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C.  HARDWARE  DESCRIPTION 

The  hardware  arrangement-  in  block  diagram  form  is  sh.own 
in  Fig.  14.  The  only  change  to  the  RTS5 D  was  the  redesign 
and  relocation  of  the  thruster  console  to  the  pilots  seat 
where  it  can.  be  more  conveniently  utilized  by  the  operator. 
The  thruster  console  uses  "linear  movement"  potentiometers 
in  place  of  previously  used  "dial"  potentiometers  in  order 
to  more  realistically  simulate  thrust  controls  for  the  SES. 
The  thruster  console  is  shown  in.  Fig.  15. 

D.  SOFTWARE  DESCRIPTION 

The  remodeled  software  package  was  a  FORTRAN  IV  digital 
program  which  consisted  of  a  main  graphics  program  and  three 
major  subroutines.  It  was  necessary  to  make  a  third  sub¬ 
routine  out  of  the  equations  of  motion  computation  that  had 
previously  resided  in  the  main  program  prior  to  making  any 
changes  to  these  equations.  It  was  discovered  when  the 
first  program  modifications  were  attempted  that  the  compiler 
for  the  XDS-9300  computer  was  operating  at  its  limit  with 
the  program  as  it  then  existed,  therefore,  before  changes  to 
the  equations  could  be  implemented,  the  program  structure 
had  to  be  modified.  The  failure  mode  was  eliminated  from 
the  new  model  because  of  initial  assumptions  that  the  thrust 
would  be  all  applied  at  the  same  point..  Thrust  failures 
should  not  be  considered  until  the  model  includes  broaching 
dynamics  which  is  a  major  source  of  thrust  loss.  The  flow 
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chart  for  the  remodeled  RTSSD  is  shown  in  Fig.  Id.  The  com¬ 
plete  multiplexing  algorithm  is  shown  in  a  flow  chart  in 
Fig.  17. 
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A  -  Thruster  T7 
B  -  Thruster  T8 
C  -  Thruster  T9 
D  -  Thruster  TlO 
E  -  Speed  Command  Control 
F  -  Run,  Hold,  Restart 

Figure  15 


Thruster  Console 


Figure  16 


RTS5D  Mod  Program  Flow  Chart 
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DESCRIPTION  OF  INFO(I) 


INFO(l)  u  DYNAMICS 
INFO(2)  u  DYNAMICS 
INFO (3)  v  DYNAMICS 
INFO (4)  v  DYNAMICS 
INFO(5)  r  DYNAMICS 
INFO(6)  r  DYNAMICS 
INFO(7)  p  DYNAMICS 
INFO(8)  p  DYNAMICS 
INFO (9)  q  DYNAMICS 
INFO (10)  q  DYNAMICS 
INFO(ll)  POSITION  MAX  VALUES 
INF0(12 )  POSITION  MIN  VALUES 


INF0(13)  VELOCITY  MAX  VALUES 
INFO (14)  VELOCITY  MIN  VALUES 
INFO (15)  ACCELERATION  MAX  VALUES 
INF0(l6)  ACCELERATION  MIN  VALUES 
INFO (17)  NAVIGATION  DATA 
INF0(18)  EFFECTOR/THRUSTER  DATA 
INF0(19)  SEA  TRACK  NAV  PLOT 
INFO (20)  SEA  TRACK  NAV  PLOT 
INF0(21)  PRESENT  POSITION  VARIABLES 
INFO (22)  PRESENT  VELOCITY  VAREABLES 
INFO (23)  PRESENT  ACCELERATION  VAR. 


Figure  17 

RTS5D  MOD  DISPLAY  MULTIPLEX  ALGORITHM 
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VI.  RTS5D  MODS  I  AND  II  RESPONSE  CHARACTERISTICS 


RTS5D  MOD  I  includes  negative  drag  effects  in  a  sea  state, 
simplified  propulsion  dynamics,  speed  control  and  a  broach 
condition  warning  flag.  RTS5D  MOD  II  is  identical  to  MOD  I 
except  that  it  includes  the  correction  factor  for  thrust 
effector  angle  developed  in  Chapter  IV,  Section  E. 

The  validation  of  the  modified  RTSSD  was  accomplished 
in  the  same  manner  as  the  RTSSD.  CSee  Chapter  1 1 1)  .  The 
results  are  displayed  in  the  following  tables.  A  complete 
560°  turn  comparison  of  the  DBSIM5D,  the  RTSSD  and  the  modi¬ 
fied  RTSSD  at  50  knots,  15°  thruster  angle  in  shown  in  Fig. 

18.  Additionally,  a  comparison  of  transient  response  was 
undertaken  and  the  time  of  first  peak  for  the  variable  v,  r, 

$  and  a  for  each  model  was  graphed  in  Figs.  19-21  for  speeds 
of  40  ,  50  and  60  knots.  Transient  response  characteristics 
are  an  important  consideration  in  the  development  of  a  real 
time  model  and  are  discussed  in  Appendix  A.  The  transient 
response  of  the  RTSSD  models  was  faster  than  that  of  the 
data  base  model  with  the  exception  of  roll  angle  (4) .  This 
impacted  on  the  validity  of  a  real  time  solution  of  the  equa¬ 
tions  of  motion  and  the  display  of  the  results. 


TABLE  IV 


RTS5D  Mod  I 

and  Mod  II 

Performance 

Test  at  *10 

Knots 

effector 

angle 

RTS5D  Mod  1 

RTS5D  Mod 

II 

* 

in  degrees 

5® 

1st  pk 
.  52 

Qss 

.45 

1st  pk 
.40 

Qss 

.38 

10° 

.  96 

.92 

.79 

.  76 

9 

in  degrees 

9 

5° 

1st  pk 
1.21 

Qss 

1.17 

1st  pk 
1.20 

Qss 

1.18 

10° 

1.20 

1.18 

1.20 

1.18 

u 

in  ft/sec 

u 

5° 

1st  pk 
62.56 

Qss 

62.64 

1st  pk 
62.58 

Qss 

67.65 

10° 

n/a 

61.37 

n/a 

61.71 

V 

in  ft/sec 

V 

5° 

1st  pk 
2.62 

Qss 

2.47 

1st  pk 
2.33 

Qss 

2.27 

10° 

3.53 

3.52 

3.23 

3.20 

r 

in  degrees/sec  r 

5° 

1st  pk 
.69 

Qss 
.  59 

1st  pk 
.  57 

Qss 

.  54 

0 

o 

t — 4 

1.18 

1.22 

1 . 04 

1.09 

62 
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TABLE  V 


RTS5D  Mod  I  and  Mod  II  Performance  Test  at  50  Knots 


effector 


angle 

RTS5D  Mod 

i 

RTS5D  Mod 

ii 

* 

in 

degrees 

* 

5° 

1st  pk 
.  50 

Qss 

.39 

1st  pk 
.43 

Qss 

.39 

10° 

.90 

.82 

.  83 

.  78 

15° 

1.31 

1.36 

1.24 

1.19 

0 

in 

degrees 

0 

5° 

1st  pk 
1.23 

Qss 

1.19 

1st  pk 
1.23 

Qss 

1 .20 

h-» 

O 

o 

1.24 

1.19 

1.23 

1.20 

15° 

1.24 

1.18 

1.23 

1.19 

u 

in 

ft/sec 

U 

5° 

1st  pk 
n/ a 

Qss 

82.73 

1st  pk 
n/ a 

Qss 

82.98 

10° 

n/ a 

78.68 

n/a 

79.55 

15° 

n/a 

68.30 

n/ a 

72.75 

V 

in 

ft/sec 

V 

5° 

1st  pk 
2.57 

Qss 

2.30 

1st  pk 
2.39 

Qss 

2.29 

H-* 

O 

o 

3.48 

3.35 

3.25 

3.25 

15° 

4.13 

4.29 

3.98 

4.00 

r 

in 

degrees/sec 

r 

5° 

1st  pk 
.  57 

Qss 

.39 

1st  pk 
.49 

Qss 

.41 

10° 

.96 

.  85 

.88 

.86 

15° 

1.31 

1.62 

1.25 

1.41 
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TABLE  VI 


RTS5D  Mod  I  and  Mod  II  Perforrr.ar.ee  Test  at  5  0  Knots 


effector 

angle  RTS5D  Mod  1  RTS  SI)  Mod  II 

■p  in  degrees  p 


S° 

S° 

5° 


1st  pk 

Qss 

1st  pk 

Qss 

.  53 

.39 

.46 

.41 

e 

in 

degrees 

9 

1st  pk 

Qss 

1st  pk 

Qss 

1.23 

1.20 

1.24 

1.20 

u 

u 

1st  pk 

Qss 

1st  pk 

Qss 

n/a 

100.21 

n/a 

99.59 

O  X 


.0 

ox  A 

A 


_  V*  0*  0? 

.  A  A 

$ 


* 


0  * 


4 

M 


•  V 


•1  O 


*  *A  •  &  O 

O  >r'  .'•  *  a 

0  ™  **  £  z  0° 


O  O 


u  =  50  knots 
z  =  15°  (effector  angle) 

10  second  intervals 

O  RTS5D 

*  DBSIM5D 

•  RTS5D  MOD  I 
X  RTS5D  MOD  II 

Figure  18 

9 

360  Turn  Comparison  of  Four  Models 
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VII.  CONCLUSIONS 


The  RTS5D  has  been  shown  in  Section  VI  to  be  a  viable 
input/output  model  of  a  3K-SES  by  closely  matching  its  output 
characteristics  for  a  given  input  to  those  of  the  data  base 
model . 

The  refinement  of  the  right  side  of  the  simplified  equa¬ 
tions  of  motion,  especially  the  inclusion  of  the  negative 
drag  characteristic  in  a  sea  state  has  made  the  model  more 
accurate  over  a  wide  speed  range  and  gives  the  operator,  upon 
program  initialisation,  the  ability  to  choose  different  drag 
characteristics  caused  by  sea  state. 

The  model  was  tested  in  a  negative  drag  region  with 
automatic  speed  control  and  simplified  propulsion  dynamics 
and  it  was  demonstrated  that  by  properly  adjusting  controller 
gains  that  the  output  characteristics  of  the  DBSIM5  9  were 
closely  matched. 

Previously  conducted  thrust  failure  analysis  with  man 
generated  corrective  action  has  been  shown  to  be  invalid 
until  broach  dynamics  are  included  into  the  model.  Previous 
tests  allowed  maximum  thrust  effector  angle  correction  regard¬ 
less  of  this  consideration.  The  modified  version  of  the  RTSS D 
does  not  include  complete  broach  dynamics  but  does  warn  the 
operator  if  he  should  exceed  a  maximum  thrust  effector  angle 
for  a  given  forward  speed. 


The  3K-SES  may  be  operated  in  the  negative  drag  slope 
speed  region  approximately  (36-58)  knots  without  the  aid  of 
a  speed  controller  but  would  require  constant  manual  modula- 
t.ion  of  the  thrusters.  Since  minimum  manning  is  one  of  the 
criteria  in  3K-SES  design,  a  speed  controller  is  necessary 
to  free  operator  attention  to  other  duties  when  operating 
in  this  speed  range. 

Another  consideration  of  the  implementation  of  the  non¬ 
linear  drag  curves  is  that  for  a  given  sea  state,  two  values 
of  minimum  drag  can  be  determined,  one  in  Region  IV  and  one 
between  Region  II  and  Region  III.  Knowledge  that  drag  can 
be  lower  in  a  particular  operating  region  for  fuel  conserva¬ 
tion  is  of  prime  importance. 

It  is  recognized  that  a  time  delay  exists  in  the  solution 
of  the  equations  of  motion  and  a  time  delay  exists  in  the 
graphical  display  of  the  solutions  to  the  operator.  Real 
time  interaction  between  the  man  in  the  loop  and  the  computer 
is  achieved  by  setting  the  integration  step  size  equal  to 
th.e  computer  loop  iteration  time.  Real  time  considerations 
are  further  discussed  in  Appendix  A. 


70 


i 


\ 


VIII .  RECOMMENDATIONS 


The  following  recommendations  are  made  concerning  future 
improvements  in  the  RTS5D  model, 

1.  Further  refinement  of  the  right  hand  side  of  the 
equations  of  motion  is  necessary  in  order  to  more  closely 
match  the  characteristics  of  the  data  base  model. 

2 .  Implement  software  broach  dynamics  in  the  model  with 
thrust  applied  differentially  at  different  lever  arms  Csee 
Fig.  3)  so  that  thrust  failures  and  effectiveness  of  correc¬ 
tive  action  may  be  evaluated. 

3.  Rewrite  the  program  such  that  the  AGT-10's  would 
provide  asyncronous  parallel  graphics  Cexisting  equipment) 
or  implement  the  model  on  a  more  modern  system  which  uses 
parallel  processing  and  improved  graphics. 

4.  Introduce  the  sixth  degree  of  freedom  Cheave)  as  a 
parallel  operation.  Since  heave  time  constant  is  much  faster 
than  the  rest  of  tire  system  dynamics,  real  time  heave  analysis 
could  be  conducted  simultaneously  if  cross  coupling  is  neg¬ 
lected. 

5.  The  system  is  designed  for  man  in  the  loop  interface, 
so  further  operator  testing  is  recommended  as  a  valuable 
source  of  information  for  future  refinement  of  the  model. 
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APPENDIX  A 


REAL  TIME  ANALYSIS 

T.  S.  Nelson,  in  Ref.  1,  accomplished  real  time  solution 
of  SES  equations  of  motion  with  the  use  of  a  1000  Hz  clock 
to  measure  computer  iteration  time  and  by  using  that  time  as 
an  integration  step  size  for  rectangular  integration  to  solve 
the  simplified  dynamic  equations.  He  concluded,  by  comparing 
the  solutions  of  these  equations  to  that  obtained  by  using  a 
Runga  Kutta  variable  integration  step  size,  that  the  computer 
iteration  loop  time  must  stay  below  100  ms  in  order  to  main¬ 
tain  at  least  3%  accuracy  of  solution.  He  further  developed 
a  multiplexing  algorithm  in  order  to  output  these  solutions 
to  an  operator  in  the  form  of  graphics  on  a  real  time  basis. 
This  model  can  be  invalidated  as  a  real  time  model  in  two 
ways.  First,  in  order  to  have  real  time,  the  integration 
step  size  used  must  be  equal  to  the  computer  iteration  time 
and  in  this  model,  because  of  the  multiplexing  algorithm, 
each  loop  time  is  different.  The  integration  step  size, 
dt  determined  by  measuring  loop  N  is  used  to  solve  the  equa¬ 
tions  in  loop  N+l,  which  has  a  different  loop  time,  there¬ 
fore  the  solutions  are  not  exactly  real  time.  Secondly, 
the  multiplexing  scheme  makes  the  model  into  a  sampled  data 
system  with  sample  time  (Tsl  approximately  2  seconds.  Even 
if  the  first  condition  did  not  apply  and  real  time  solutions 
to  the  equations  could  be  assumed,  then  the  output  to  the 
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operator  would  be  sampled  real  time  data.  Given  the  fastest 
time  constant  for  the  RTS5D  simplified  equations  (approxi¬ 
mately  1.3  seconds  for  pitch  rotational  motion),  it  can  be 
seen  that  the  Nyquist  Criteria  is  clearly  violated.  Yet 
even  with  these  obvious  limitations,  it  has  been  demonstrated 
that  the  RTS5D  can  closely  reproduce  the  output  characteris¬ 
tics  of  the  DBSIM5D  and  therefore  is  a  useful  model.  Surge 
speed  as  well  as  sway,  yaw  and  roll  (Figs.  21-73)  have  much 
longer  effective  time  constants.  It  is  recognized  that 
there  is  a  delay  between  the  real  time  of  any  event  and  the 
observation  of  that  event,  therefore  "Real  Time"  as  used  in 
this  report  is  a  matter  of  degree,  or  more  specifically 
how  much  delay  time  between  solution  and  display  can  be 
tolerated.  Given  the  scale  of  graphics  presented  to  the 
operator  in  this  model  and  the  speed  with  which  the  SES 
responds,  the  RTS5 D  is  an  acceptable  model,  however  it 

would  be  desirable  to  reduce  the  samDle  time  T  and  make  the 

s 

integration  step  size  equal  to  the  computer  loop  iteration 
time  for  each  iteration. 

An  algorithm  to  set  the  integration  step  size  equal  to 
the  computer  iteration  time  of  the  current  loop  was  imple¬ 
mented  as  a  software  design  change  to  the  RTS5  D  model.  This 
was  accomplished  by  defining  the  integration  step  size 
as  an  array  DELTC  N)  with  N  equaL  to  th  e  total  number  of  itera¬ 
tions  necessary  to  display  all  of  the  graphical  information 
to  the  operator.  Each  loop  iteration  time  was  Treasured  and 


stored  in.  this  array  during  the  first  pass  through,  the 
multiplexing  algorithm.  DELTCN7  was  initialized  at  taro 
so  that  time  would  not  increment  and  integration  would  not 
be  performed  during  the  first  pass  through  the  multiplexing 
algorithm.  On  subsequent  passes  through  the  multiplexing 
algorithm  the  integration  step  size  was  equal  to  the  loop 
time  in  which  it  was  used. 

this  timing  algorithm  when  tested  in  the  RTS5 D  model 
was  found  to  accurately  represent  real  time  but  the  devia- 
t ion  from  previously  used  integration  step  size  values 
caused  instability  in  both  the  pitch  and  roll  dynamics. 

This  needs  further  investigation. 

Sturgeon,  in  Ref.  8,  demonstrated  the  value  of  using 
parallel  array  processors  to  reduce  computational  time  in 
real  time  simulations.  Such  a  model  could  be  designed  for 
the  RTS5D  using  the  existing  equipment.  The  AGT-lO's  which 
are.  used  for  displaying  the  output  can  interface  with  the 
XDS-33G0  main  computer  in  two  ways.  The  multiplexing  scheme 
of  the  RTS5D  uses  a  "hand  shaking"  method  with  a  GRAPHO  or 
TEXTO  subroutine  call  in  each  loop.  It  is  precisely  the 
inclusion  of  graphics  in  each  loop  that  causes  problems  for 
real  time  presentation.  The  solution  of  the  equations  requires 
only  approximately  35  ms.  Each  GRAPHO  or  TEXTO  call,  with 
its  associated  "hand  shaking"  time,  makes  up  the  rest  of  the 
loop  time,  approximately  40  ms  more.  It  can  be  seen  that  by 
introducing  two  such  calls  in  one  loop  that  the  100  ms  time 
restriction  would  be  exceeded. 
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An  improvement  to  this  method  would  be  to  use  the  second 


interface  between  the  AGT-lQ?s  and  the  XDS-330Q  which  is 
that  the  AGT-lO's  have  a  direct  connection  to  some  memory 
locations  in  XDS-3200.  By  solving  the  dynamics  for  the  SES 
and  storing  the  data  in  these  memory  addresses  in  the  form 
of  arrays ,  the  AGT-TJrs  could  be  made  to  operate  asyncronously 
in  parallel  as  array  processors  by  using  a  process  known  as 
"cycle  stealing"  to  fetch  the  data  from  the  XDS-9300  and  to 
display  this  on  the  screen  for  the  operator.  This  process 
would  have  the  advantage  of  reducing  and  making  more  uniform 
the  computer  iteration  time  used  as  the  step  size  for  the 
rectangular  integration  in  the  solution  of  the  SES  equations 
of  motion.  The  sample  time  would  also  be  reduced  because  all 
of  the  time  wasted  in  the  "hand  shaking"  is  eliminated.  It's 
drawback  is  that  the  output  is  not  trueiy  real  time  because 
of  the  asyncronous  operation  of  the  AGT-10Ts.  Kowever,  the 
reduced  sampled  data  time  and  computer  iteration  time  would 
make  it  a  closer  approximation  of  real  time  than  presently 
exists , 
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All 
A21 
A22 
A  31 
A?3 
A 

AM 

ABSYD 

A23YCD 

AD 

A1CCA 

AM 

APITCH 

A PRINT 

ARR 

A  RATE 

ARCLL 

AST  CP 

AXST 

AYST 

BETA 

SCAT 


APPENDIX  E 

RTS5D  MODIFIED  PROGRAM  NOMENCLATURE 

A/D  trunk  500  voltage 
A/D  trunk  $0 1  voltage 

Added  mass  coefficient  in  yaw  moment  equation 
A/D  trunk  502  voltage 

Added  mass  coefficient  in  roll  moment  equation 
Added  mass  coefficient  in  pitch  moment  equation 
A/D  trunk  503  voltage 

y 

v 

‘0 

array  of  A/D  lines 
speed  line  index  counter 
mass 

pitch  angle  in  degrees 

f  iterations  between  print  execute  commands 

Real  IARR  (l) 

turn  rate  in  degrees 

roll  angle 

performance  index  limit 
initial  condition  X-coordinate  on  sea  track 
initial  condition  Y-coordinate  on  sea  track 
drift  angle 

subroutine  to  generate  perspective  of  3SS 
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CA33T 

CDX 
CDY 
CDZP 


(i) 


DAL 
DELT 
DELTA 
DEC  CM 
DEC CM  1 


DEC  CM  2 


DTIME 

DTIMPLT 

DUG  CM 

DXPLCT 

DYPLCT 

DYREPET 

EFLG 

ETIME 

EYRE PET 

EXPLOT 

EYPLOT 

FI 

F2 

F3 


automatic  speed  control  mode  indicator 

cow  real  forces  lumped  coefficient 

drag  forces  lumped  coefficients 

sway  forces  lumped  coefficient 

sidewall  rolling  moment  lumped  ceofficient 

digital  to  analog  call 

iteration  loop  time 

commanded  iteration  loop  time 

total  change  in  fuel  command  from  speed  comtroller 
change  in  fuel  command  due  to  proportional  speed 
controller 

change  in  fuel  command  due  to  proportional  speed 
controller 
v  time 

scaled  value  for  v 
difference  between  ucom  and  u 
same  as  DTIMPLT 
scaled  value  for  v 
restart  value  for  v 
error  flag  for  TEXTC  call 
r  time 

restart  value  for  r 
scaled  r  time 
scaled  r 

stem  buoyancy  force 
bow  seal  pitch  force 
bow  buoyance  force 
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fc  cy 

9 1  c  Onirnan  d 

FSP 

port  sidewall  buoyancy  force 

FS3 

starboard  sidewall  buoyancy  force 

■71/n  «  rjT* 

v  time 

FTII-TLT 

• 

scaled  v  time 

FXPLCT 

same  as  FTIITLT 

FYPLCT 

• 

scaled  v 

PYRE PET 

restart  value  for  v 

n 

vX 

gravitational  acceleration 

C-RAPHO 

subroutine  to  project  image  (non  alpha-nuneri 

GTH'.'E 

R  time 

GTIMPLT 

♦ 

scaled  R  time 

GXFLOT 

same  asGTITTLT 

GYPLCT 

scaled  R 

GYRE PST 

restart  R  value 

H(i) 

dynamics  maximum  and  minimum  array 

HEAD 

craft  heading 

HOLD 

subroutine  to  freeze  display/progran 

HTIHE 

u  time 

HTISfiPLT 

scaled  u  time 

HYREPET 

analog  restart  y  coordinate  on  sea  track 

HY3ESX 

analog  y  position  of  craft 

HYSESY 

analog  x  position  of  craft 

IARR(l) 

digitized  mode  selector 

IARR(2 ) 

digitized  effector  angle 

IARR(3) 

digitized  T7 

78 


IARR(U) 

:arr(5) 

:arr(6) 

IARR(7) 
ICO 
ICOA 
I CRAFT 

TTTO 

1 1  PR 
1 1  PRD 

I  IFRS 
IIPRF 
IIPRG 

I I  ?RH 
IIPRC 
IIFRP 
IIPRQ 
IIPRR 
I  IFRS 
IIX 
IJ 
IJA 

I JAD 
I  JAR 
I  JAP 
I  JAG 
UAH 


digitized  T9 
digitized  T9 
digitized  TlO 

digitized  surge  velocity  command 

multiplex  graphics  index 

speed  line  index 

vehicle  designator 

error  flag 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

time  history  limit 

print  counter 

craft  perspective  counter 

craft  perspective  counter 

time  history  counter 

time  history  counter 

time  history  counter 

time  history  counter 

time  history  counter 
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IJAO 


IJAO 

time 

history 

counter 

I  JAP 

time 

history 

counter 

IJAQ 

time 

history 

counter 

I  JAR 

time 

history 

counter 

I JAS 

time 

history 

counter 

IJD 

time 

history 

counter 

IJE 

time 

history 

counter 

IJF 

time 

history 

counter 

IJG 

time 

history 

counter 

IJH 

time 

history 

counter 

IJO 

time 

history 

counter 

IJP 

time 

history 

counter 

IJQ 

time 

history 

counter 

IJR 

time 

history 

counter 

IJS 

time 

history 

counter 

IL 

graphics 

console 

display 

array 

I  LAI 

graphics 

console 

display 

array 

ILA2 

graphics 

console 

display 

array 

ILA3 

graphics 

console 

display 

array 

ILA4 

graphics 

console 

display 

array 

ILA5 

graphics 

console 

display 

array 

ILA6 

graphics 

console 

display 

array 

ILA7 

graphics 

console 

display 

array 

ILA8 

graphics 

console 

display 

array 

ILA9 

graphics 

console 

display 

array 

ILA10 

graphics 

console 

display 

array 
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ILA11 

graphics 

console 

display 

array 

I  LAI  2 

graphics 

console 

display 

array 

I  LAI  3 

graphics 

console 

display 

array 

I  LAI  4 

graphics 

console 

display 

array 

I  LAIS. 

graphics 

console 

display 

array 

I  LAIS 

graphics 

console 

display 

array 

I  LAI  7 

graphics 

console 

display 

array 

I LAI  3 

graphics 

console 

display 

array 

I  LA  13 

graphics 

console 

display 

array 

ILA23 

graphics 

console 

display 

array 

ILA21 

graphics 

console 

display 

array 

I  LA?  2 

graphics 

console 

display 

array 

ILA23 

graphics 

console 

display 

array 

ILA2  4 

graphics 

console 

cispl  ay 

array 

ILA2S 

graphics 

console 

display 

array 

ILA25 

graphi cs 

console 

display 

array 

ILA2  7 

graphics 

con  sole 

display 

array 

ILA2  3 

graphics 

con  sole 

cispl ay 

array 

ILA.29 

graphics 

console 

display 

array 

LA? 

graphics 

console 

displ ay 

array 

I  PACK 

subroutine  to  load  displa 

arr< 

I  PLOT 

graphi cs 

console 

display 

array 

IPR 

time  history  counter  limit 

I  PRD 

time  hist 

ory  coun 

ter  limi 

t 

IPRE 

IPRF 

IPRG 

IPRH 

IPRO 

IPRP 

IPRQ 

IPRR 

IPRS 

ISWA 

I  TEXT 

IVEIWA 

IVIEWAA 

IVIEWB 

IVIEWC 

IVIEWD 

IVIEWE 

IVIEWF 

IVIEWG 

IVIEWH 

IVIEWO 

IVIEWP 

IVIEWQ 

IVIEWR 

IVIEWS 

IVIEWZ 


time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
time  history  counter  limit 
display  array 

subroutine  to  display  alpha-numeric  data 

craft  perspective/north  road  array 

sea  track  array 

sea  track  axis  array 

u  time  history  array 

v  time  history  array 

r  time  history  array 

v  time  history  array 

r  time  history  array 

u  time  history  array 

p  time  history  array 

q  time  history  array 

p  time  history  array 

q  time  h.i  story  array 

analog  sea  track  array 

time  history  axis  array 
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J 

JD 

U  — 

jp 

JG 

JH 

JC 

JP 

JQ 

JH 

JS 


L3 

LD 

LP 

12 

*’< 

NAD 

ONET'CT 

CO 

CTI?"E 
CTI'TP'’  T 
CXPLCT 
CYPLCf 
CYRSPET 


x  axis  moment  of  inertia 

fine  history  counter 

tine  history  counter 

tine  history  counter 

time  histroy  counter 

time  history  counter 

time  history  counter 

tine  history  counter 

tine  history  counter 

tine  history  counter 

tine  history  counter 

tine  history  counter 

proportional  speed  controller  gain 

intregal  speed  controller  gain 

lever  arm  for  buoyancy  force 

draft  of  sidewall 

side  thrust  lever  arm 

lever  arm  for  sway  drag  force 

clock  interrupt  count 

AD  trunk  line  array 

subroutine  to  calculate  SSS  dynamics 

lever  arm  for  thrust  side  component 

• 

p  time 

scaled  p  tine 
sane  as  CTINPLT 
scaled  p 
restart  p  value 


33 


roll  rate 


?5 

plenum  pressure 

PDCT 

roll  acceleration 

?dct:-: 

max  roll  accelera' 

PHI 

roll  angle 

PKAX 

max  roll  rate 

PSI 

craft  heading 

PTIHE 

q  time 

p  (Tl  T  '  ro  T  T1 

• 

scaled  q  time 

FXPICT 

same  as  PTI'TLT 

• 

scaled  q 

FYR2PST 

restart  q  value 

'4 

pitch  rate 

;dct 

pitch  acceleratioi 

3DCT!.: 

max  a 

QMAX 

max  q 

QTir.z 

p  time 

tit  * 'rC T  rn 

'ill  .-JT L.J. 

scaled  p  time 

QXrICT 

same  as  .JTIHPLT 

3YPLCT 

scaled  p  value 

QYREFET 

restart  p  value 

R 

craft  turn  rate  r 

P.DCT 

RDCTH 

max  r 

RZADCLCCK  subroutine  to  sample  real  time  clock 
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RHC 

rhax 

j  - -n? 

RTII.TLT 

RX 

RXPLCT 


subroutine  to  reset  analog  computer 
density  of  water 
maxinum  yaw  rate  of  craft 
q  rate 

scaled  q  time 
IX 

same  as  RTITIPLT 


RY 


I 


V 


RYPLCT 

RYREPET 

RZ 

SEAST 

3ESX 

3E3Y 

SF1 

3F14 

SLIP 

SPDLin 

SPEED 

3TARTCLCCX 

SrpT\*TT 

X  a*  *■  -j—i 

STIXPLT 
SXPLOT 
SYPLCT 
3 X RE PET 


scaled  q  value 
restart  q  value 

T 

x  n 
u 

sea  state 
scaled  X 
scaled 

o 

scale  factor  in  time  history  arrays 

scale  factor  in  time  history  arrays 

Beta  influence  coefficient 

speedline  lower  screen  limit 

Craft's  velocity  in  knots 

subroutine  to  start  clock 

analog  time 

scaled  analog  time 

HYSESX 

HYSE3Y 

restart  value  for  HYSESX 


3YRSPET 


restart  value  for  HYSESY 


otal  thrust 


T10 

thruster  nuhcer  4 

TIOHAX 

max  T10  value 

T? 

thruster  number  1 

T7MAS 

max  T7  value 

T8 

thruster  number  2 

TSXAX 

max  value  of  T8 

T9 

thruster  number  3 

Tormx 

max  value  of  T9 

TEXTC 

subroutine  to  display  alpha-r.umer 

TDCT 

T 

THETA 

pitch  angle 

TI" 

u  time 

TIT77 

time 

TIMPLT 

scaled  TIM 

TINT 

real  clock  interrupt  frequency 

TITLE 

program  name  display  array 

TITLED 

program  name  display  array 

TITLE1 

program  name  display  array 

T»i 

manual  thrust  command 

TLAX 

max  T 

TSIDE 

sum  of  side  thrust  components 

TYA'/J 

sum  of  yaw  moment 

TT 

V 

forward  velocity  (surge) 

UCCH 

forward  velocity  (surge)  command 
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l 


V 


UDCT 


u 


UDCT'.:  max  u 

UMAX  max  u 

UUUL  scale  factor  in  speedline  generation 

V  side  velocity  (sway) 

VCD  thruster  console  pot  array 

VDOT  v 

VDCTM  max  v 

WAX  max  v 

VS  total  velocity 

'■i  heave  velocity 

'VDCT  heave  acceleration 

'.VS  width  of  how  seal 

VJRITECLCCK  subroutine  to  assign  clock  interrupts 
variable 

WW  lever  arm  for  sway  drag  forces 

X  Y 

o 

XC  X 

o 

XCDCT  X 

o 

XDX  boat  perspective  scale  factor 

X3X  boat  perspective  scale  factor 

XPLCT  display  array 

X RE PST  reset  value  for  SESX 

XST  scale  factor  for  analog  plot 

VC  V 

xo 

YCDCT  Y 

o 


37 


-HI  3H 
Y?L  CT 
YRZFET 
■'ST 
Z 

12 
Z3  ' 
zu 

Z  5 
ZA3 
Z33C 
21 

77  7-* 

O  O  £4 


point  D  vertical  position  constant 
display  array 
reset  value  for  SZSY 
scale  factor  for  analog  plot 
operational  effector  angle 
upper  effector  angle  limitation 
lower  effector  angle  limitation 
zero  angle  upper  limit 
zero  angle  lower  limit 
scale  factor  for  Z10 
effector  broach  limit  angle 
modified  Z 

fixed  step  effector  angle  input 


APPENDIX  C 

RTS5D  MODIFIED  COMPUTER  PROGRAM  LISTING 
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